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ABSTRACT 

/. 

A  new  Ildar  technique  is  proposed  for  determining  aerosol 
optical  parameters.  It  consists  in  recording  the  lidar  returns  at 
three  fields  of  view  to  measure  the  differences  caused  by  multiple 
scatterings.  The  inversion  algorithm  uses  that  information  to  obtain 
the  aerosol  forward-scattering  coefficient  through  simple  algebraic 
formulas.  Sample  solutions  are  calculated  for  simulated  lidar  signals 
and  compared  with  the  results  derived  from  the  standard  single¬ 
scattering  lidar  equation.  The  agreement  with  the  true  values  is  very 
good  and  performance  is  consistently  better  than,  or  at  least  as  good 
as,  that  of  the  single-scattering  method.  The  proposed  technique 
requires  neither  a  boundary  value  nor  a  backscatter-to-extinction 
relation,  which  makes  such  lidar  measurements  completely  self- 
sufficient  for  the  remote  determination  of  the  range-resolved  forward¬ 
scattering  coefficient  of  atmospheric  aerosols.  ^ 

RESUME 


Une  nouvelle  technique  lidar  pour  la  determination  de  param3tres 
optiques  des  aerosols  est  proposee.  Elle  consiste  3  enregistrer  les  / 

retrodif fusions  lidar  sous  trois  champs  de  vision  angulaire  pour  mesu- 
rer  les  differences  provenant  des  diffusions  multiples.  L'algorithme 
d' inversion  utilise  ces  donnees  pour  obtenir  par  des  operations  alge- 
briques  simples  le  coefficient  de  diffusion  vers  l'avant  des  aerosols. 

Des  exemples  de  calcul  sont  effectues  pour  des  signaux  lidar  simuies  et 
compares  aux  resultats  obtenus  3  partir  de  liquation  du  lidar  en  dif¬ 
fusion  simple.  Le  recoupement  avec  les  valeurs  exactes  est  tres  bon  et 
les  performances  sont  meilleures  ou  du  moins  aussi  bonnes  que  celles  de 
la  methode  de  diffusion  simple.  La  methode  suggeree  requiert  ni  valeur 
aux  limites  ni  relation  entre  les  coefficients  d* extinction  et  de 
retrodif fusion,  ce  qui  la  rend  tout  £  fait  autosuf f isante  pour  dete 
ner  3  distance  la  distribution  spatiale  du  coefficient  de  diffusion 
vers  l'avant  des  aerosols  atmospheriques . 
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1.0  INTRODUCTION 


A  major  characteristic  of  modern  defence  systems  is  the  use  of 
electro-optic  (E0)  devices  for  search,  detection,  tracking,  guidance 
and  control.  The  performance  of  these  systems  Is  very  good  in  clear 
visibility  but  it  can  be  severely  degraded  by  poor  weather  conditions, 
especially  by  the  presence  of  aerosol  particles.  The  principal  diffi¬ 
culty  in  the  tactical  deployment  of  these  systems  is  the  great  vari¬ 
ability  of  aerosol  composition  and  concentration,  which  makes  perfor¬ 
mance  assessment  at  any  time  in  any  weather  very  unreliable.  There¬ 
fore,  there  is  a  need  for  an  instrument  that  could  remotely  measure  in 
real  time  and  in  situ  the  aerosol  extinction  coefficient  in  the  imme¬ 
diate  environment  where  the  E0  systems  are  to  be  deployed. 

Lidar  appears  as  an  attractive  solution  to  this  problem  of 
determining  the  aerosol  extinction  coefficient  In  an  inhomogeneous 
medium.  However,  the  technique  still  involves  a  number  of  difficul¬ 
ties.  The  standard  single-scattering  lidar  theory  assumes  that  the 
detected  photons  have  encountered  only  one  scattering  before  being 
collected.  Moreover,  the  solution  method  requires  a  relation  between 
the  backscatter  and  extinction  coefficients,  and  the  specification  of 
boundary  value  of  the  extinction  coefficient,  preferably  at  the  outer 
end  of  the  lidar  range.  These  data  cannot  in  general  be  derived  from 
the  measured  lidar  signals;  most  often,  they  are  guessed  or  estimated 
through  very  limited  Independent  information.  Because  of  these  prob¬ 
lems,  the  single-scattering  lidar  solutions  are  subject  to  important 
errors  and  uncertainties  (Refs.  1-5)  that  can  make  tactical  applicatio 
of  lidar  unfeasible  under  general  weather  conditions. 
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This  report  proposes  a  Ildar  technique  that  uses  the  information 
carried  by  the  multiple-scattering  contributions  to  relax  the  limita¬ 
tions  imposed  by  the  6ingle-scattering  equation.  The  method  entails 
the  simultaneous  recording  of  lidar  returns  at  three  fields  of  view  but 
it  requires  neither  a  boundary  value  nor  a  backscatter-to-extinction 
relation.  The  inversion  algorithm  is  simple,  it  is  similar  to  the 
well-known  slope  method  but  applied  to  the  ratios  of  the  lidar  signals 
at  the  different  fields  of  view.  Sample  calculations  are  performed 
with  simulated  signals  and  the  results  are  promising.  In  a  recent 
article,  Sassen  and  Petrilla  (Ref.  6)  hinted  that  multiscattering  con¬ 
tributions  can  indeed  contain  recoverable  information. 

Chapter  2.0  summarizes  the  method  of  generation  of  the  simulated 
lidar  signals.  The  resulting  single-scattering  lidar  solutions  are 
given  and  discussed  in  Chapter  3.0.  Finally,  Chapter  4.0  describes  the 
proposed  multiscattering  lidar  concept  and  discusses  the  solutions 
obtained. 

This  work  was  performed  at  DREV  between  January  and  May  1986 
under  PCN  33B15,  Lidar  Determination  of  Atmospheric  Parameters. 

2.0  LIDAR  SIMULATION 


2.1  Model 


In  analyzing  the  lidar  solutions  for  the  aerosol  extinction 
coefficient,  it  is  important  to  compare  the  results  with  the  true 
coefficient.  This  is  difficult  to  achieve  experimentally  since  the 
lidar  gives  a  range-resolved  distribution  of  the  extinction  coefficient 
which  contains  information  on  a  finer  scale  than  generally  available 
from  independent  simultaneous  measurements.  Transmittance  data  are 
useful  but  they  give  only  the  range-integrated  extinction  coefficient. 


•  ■«***,  j-  -  v y  ^ 
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Vertical  profiles  of  visibility  and  particle  size  distributions  have 
been  obtained  during  lidar  experiments  at  the  same  site  (Refs.  7  and 
8).  The  visibility  data  can  be  transformed  to  extinction  coefficients 
but  the  wavelength  is  different  than  that  of  the  lidar.  Similarly, 
extinction  profiles  can  be  calculated  from  the  size  distributions  but 
the  results  are  subject  to  unknown  errors  because  of  instrumental  range 
limitations  on  particle  sizes  and  uncertainties  in  the  complex  refrac¬ 
tive  index  of  the  particles.  Although  such  measurements  are  indis¬ 
pensable  and  constitute  the  ultimate  test,  they  cannot  provide  detailed 
comparisons  under  varied  conditions.  A  convenient  alternative,  suit¬ 
able  for  comparing  solution  methods,  is  to  numerically  generate  lidar 
signals  for  given  profiles  of  extinction  and  backscatter  coefficients. 
In  most  studies  published  to  date  (e.g.  Refs.  1-4) ,  the  model  used  to 
compute  the  lidar  signals  is  the  single-scattering  lidar  equation 
itself.  Hence,  such  practical  effects  as  the  dependence  on  the  receiv¬ 
er  field  of  view  and  multiscattering  contributions  cannot  be  analyzed. 

For  the  present  application,  the  lidar  signals  are  generated 
through  solutions  of  a  multiscattering  propagation  model  developed  at 
DREV  (Refs.  9  and  10).  The  model  predicts  the  radiant  intensities  of 
the  transmitted,  forward-  and  backscattered  radiation  from  a  coherent 
laser  beam  directed  into  a  particulate  medium.  The  solutions  can  be 
calculated  for  Inhomogeneous  clouds  and  they  take  into  account  multiple 
scatterings.  Validation  was  performed  against  laboratory  data  obtained 
at  three  wavelengths  (0.63,  1.06  and  10.6  pm)  in  water  droplet  clouds 
(Ref.  10).  The  agreement  is  very  good,  especially  for  the  backscatter 
solutions  which  are  of  interest  here.  The  predictions  also  agree 
within  5Z  with  Monte  Carlo  backscatter  calculations  carried  out  at 
1.06,  3.5  and  10.5  pm  for  various  generic  aerosol  models  (Ref.  10). 
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2. 2  Cloud  Specifications 


The  aerosol  clouds  are  modeled  from  the  Air  Force  Geophysics 
Laboratory  (AFGL)  generic  rural,  urban  and  maritime  aerosols  (Ref.  11). 
These  aerosols  are  representative  of  atmospheric  conditions.  For  the 
present  application,  they  are  mixed  into  various  configurations  to 
create  a  wide  range  of  inhomogeneities.  All  clouds  are  assumed  to  be 
1  km  in  depth. 

The  first  cloud  studied  has  a  homogeneous  composition,  the 

maritime  aerosol  model  at  99 %  relative  humidity  was  chosen.  This  gives 

a  ratio  of  the  backscatter  to  extinction  coefficients  8/a  of  0.038  at 

e 

the  1.06-pm  wavelength  of  the  assumed  lidar.  The  density  is  varied  to 
give  a  saw  tooth  extinction  profile,  the  extinction  coefficient  ag  is 
1  km-1  at  the  cloud  boundaries  and  6  km-1  at  the  center.  The  cloud 
optical  depth  is  3.5.  Figure  1  shows  the  and  p/ag  profiles  for  this 
case. 


The  second  cloud  is  made  more  complex  by  varying  the  aerosol 
composition  to  modulate  the  p/ag  ratio.  The  spatial  variation  of  this 
ratio  constitutes  an  important  difficulty  for  lidar  inversion.  Urban 
and  maritime  aerosols  are  mixed  in  decreasing  and  Increasing  proportion 
respectively  over  the  first  half  of  the  cloud,  and  maritime  and  rural 
aerosols  are  similarly  mixed  over  the  second  half.  The  particles  are 
100%  urban  at  the  near-end  boundary,  100%  maritime  at  the  center  and 
100%  rural  at  the  far-end  boundary.  The  relative  humidity  is  70% 
everywhere.  This  produces  a  p/ag  at  1.06  pm  that  linearly  increases 
from  0.0096  to  0.0396  and  then  decreases  to  0.0210.  Thus,  the  p/ag 
ratio  changes  smoothly  by  a  factor  of  approximately  4  in  the  first 
half,  and  2  in  the  second  half.  The  extinction  profile  is  chosen  the 
same  as  in  Fig.  1.  Figure  2  illustrates  the  ag  and  p/ag  distributions 
for  this  case. 
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RANGE  IN  CLOUD(km) 

FIGURE  1  -  Profiles  of  the  extinction  coefficient  a  and  backscatter- 

e 

to-extinction  ratio  p/afi  at  a  wavelength  of  1.06  pm  for  a 
1-km  deep  cloud  consisting  of  the  generic  maritime  aerosol 
model  (Ref.  11)  at  99%  relative  humidity. 


00  02  04  06  08  10 


RANGE  IN  CLOUD(km) 

FIGURE  2  -  Same  as  Fig.  1  but  for  a  cloud  consisting  of  a  mixture  of 
the  generic  urban,  maritime  and  rural  aerosol  models 
(Ref.  11)  at  70%  relative  humidity. 
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A  more  complex  configuration  is  studied  by  Juxtaposing  regions 
of  different  aerosol  composition.  Figure  3  illustrates  the  structure. 
The  three  regions  have  each  a  saw  tooth  extinction  profile  but  they  are 
made  respectively  of  urban,  maritime  and  rural  aerosols  at  702  relative 
humidity.  The  p/a  ratio  thus  varies  spatially  in  sharp  steps  from 
0.0096,  up  to  0.0396  and  down  to  0.0210.  The  expected  difficulty  here 
is  the  more  complex  structure  and  the  abrupt  changes  in  the  p/ag 
profile.  Another  example  of  similar  complexity  is  shown  in  Fig.  4 
where  rural  and  maritime  aerosols  are  made  to  alternate  with  varying 
concentrations. 

The  clouds  depicted  in  Figs.  1-4  were  not  chosen  to  model  spe¬ 
cific  or  realistic  conditions  but  to  provide  a  range  of  variations  for 
both  the  extinction  coefficient  and  the  backscatter-to-extinction  ratio 
that  can  be  representative  of  natural  aerosols. 

The  lidar  returns  from  these  clouds  are  computed  with  the  propa¬ 
gation  model  outlined  In  the  preceding  section.  It  is  assumed  that  the 
transmitter  and  receiver  share  the  same  optical  axis  and  are  located  at 
the  same  distance  from  the  cloud.  The  source  is  a  nominal  1-W,  1,06-pm 
wavelength  and  100-mm  diameter  collimated  laser  beam.  The  receiver  has 
an  aperture  of  250  mm  and  an  adjustable  field  of  view.  The  range- 
resolved  lidar  backscattered  power  P(R)  is  calculated  per  unit  length 
so  that  the  power  received  from  range  R  for  a  pulsed  laser  of  average 
pulse  power  Pq  In  watts  and  pulse  duration  x  is  given  by  PqP(R)cx/2, 
where  c  is  the  speed  of  light  and  cx  is  expressed  in  km.  For  example. 
Fig.  5  shows  the  lidar  signals  P(R)  obtained  from  the  cloud  defined  in 
Fig.  3  for  a  cloud-to-receiver  distance  of  1  km  and  fields  of  view 
(half  angle)  of  1,  20,  50  and  100  mrad.  The  lidar  returns  increase 
with  field  of  view  because  of  multiscattering  contributions. 
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FIGURE  5  -  Example  of  calculated  Ildar  signals  at  receiver  fields  of 
view  (half  angle)  of  1,  20,  50  and  100  mrad  for  the  cloud 
defined  In  Fig.  3  and  for  a  cloud-to-receiver  distance  of 
1  km.  The  source  Is  a  1-W,  1.06-pm  and  100-mm  diameter 
collimated  laser  beam  and  the  receiver  aperture  diameter  is 
250  mm.  The  calculations  were  performed  with  the  propaga¬ 
tion  model  of  Ref.  10. 


3.0  SINGLE-SCATTERING  LIDAR  INVERSION 


3.1  Solution  Methods 


Although  the  object  of  this  report  is  to  introduce  a  new  Ildar 
Inversion  algorithm  based  on  recovering  the  information  contained  in 
the  multiscattering  contributions,  It  is  important  to  compare  results 
with  solutions  derived  from  the  well-known  single-scattering  lidar 
equation.  The  single-scattering  lidar  solution  methods  are  reviewed 


briefly  and  the  calculation  results  for  the  clouds  defined  in  Figs.  1-4 
are  presented  and  discussed. 


-\  -\ 
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According  to  the  single-scattering  theory,  the  range-resolved 
lidar  backscattered  power  per  unit  length  is  given  by 


ps(R)  “  p0  —  P(R)  exp  {-  2 


R 


/  ae(r)dr } , 


[1] 


where  Pq  is  the  transmitter  power,  A  is  the  area  of  the  receiver 

aperture,  R  is  the  range,  p  is  the  backscatter  coefficient,  and  a  is 

6 

the  extinction  coefficient.  A  more  convenient  signal  variable  is  the 
logarithmic  range-adjusted  power  defined  as 


S(R) 


In 


R2P  (R)' 
s 


"r2 
.  R 


2  P 
m  m 


[2] 


where  R  and  P  are  unspecified  normalization  constants.  In  terms  of 
ram 

S(R),  eq.  1  can  be  rewritten  in  the  following  differential  form: 


dS  _  1  dp 
"3r  *  ]T  cRT 


2  a  • 
e 


[3] 


Equation  3  relates  one  measured  function  S(R)  to  two  unknown  functions 
p(R)  and  ae(R).  To  solve  ae>  a  functional  relationship  must  be  given 
or  assumed  between  p  and  ag.  The  most  common  assumption  is  to  postu¬ 
late  that 


P(R)  *  C(R)  <Xg(R), 


[*] 


where  k  Is  constant.  The  function  C(R)  is  unknown  a  priori. 


The  solution  of  the  first-order  differential  eq.  3  necessitates 
the  specification  of  one  boundary  value.  To  keep  the  solution  as  gen¬ 
eral  as  possible,  we  assume  that  the  extinction  coefficient  is  known  or 
given  at  some  arbitrary  range  R^ ,  i.e.  ■  <x  (R^).  After  substitution 
of  eq.  4  for  P(R),  the  solution  of  eq.  3  (Ref.  2)  is  then  given  by 
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Equation  5  contains  the  function  C(R).  In  general,  we  cannot 
expect  to  know  C(R)  until  afi(R)  is  determined.  C(R)  is,  in  fact,  part 
of  the  solution  and  ae(R)  cannot  be  solved  for  unless  independent 
information,  in  addition  to  S(R),  is  given.  This  constitutes  one  major 
difficulty  for  the  solution  of  the  single-scattering  lidar  equation. 

One  convenient  simplification  consists  in  assuming  that  C(R)  ■  cst.  As 
shown  by  eq.  5,  the  solution  for  ag(R)  then  becomes  Independent  of  the 
value  of  C.  This  approximation  is  truly  valid  only  if  the  composition, 
geometry  and  size  distribution  of  the  aerosol  particles  are  invariant. 

It  is  emphasized  that  the  boundary  value  cannot  be  arbitrary. 
The  apparent  one-constant  arbitrariness  associated  with  the  first-order 
derivative  of  eq.  3  is  not  a  property  of  the  physical  problem  since  the 
differentiation  was  carried  out  for  convenience  only.  The  physical 
process  measured  by  the  function  S(R)  is  actually  modeled  by  eq.  1 
which  involves  no  derivatives.  Hence,  the  boundary  value  must  be 
consistent  with  S(R). 

Originally,  eq.  5  was  used  with  the  boundary  value  specified  at 
the  near-end  boundary  of  the  cloud,  i.e.  at  R^  <  R.  The  obvious  advan¬ 
tage  is  that  the  extinction  coefficient  a ^  is  easier  to  measure  or 
estimate  at  close  range  than  at  far  range.  However,  calculation 
results  (Ref.  1)  and  a  theoretical  analysis  (Ref.  5)  have  shown  that 
the  solution  given  by  eq.  5  is  unstable  If  R^  <  R.  Small  errors  cause 
the  solution  to  vanish  or  grow  indefinitely.  By  contrast,  the  same 
analyses  (Refs.  1  and  5)  have  demonstrated  that  the  solution  of  eq.  5 
is  stable  if  is  specified  at  the  far-end  boundary  or  backside  of  the 
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cloud,  l.e.  R^  >  R.  The  problem,  however,  is  to  obtain  a  reasonable 
estimate  of  o^.  This  constitutes  a  second  major  difficulty  affecting 
the  inversion  of  the  single-scattering  lidar  equation. 

3.2  Results 


To  illustrate  the  stability  problem,  we  have  applied  the  solu¬ 
tion  of  eq.  5  to  the  lidar  signal  generated  for  the  cloud  defined  in 
Fig.  1.  In  this  report,  we  make  the  parameter  k  of  eqs.  4  and  5  equal 
to  unity.  Since  C(R)  is  allowed  to  vary,  this  condition  has  no  inci¬ 
dence  on  the  generality  of  the  solutions.  For  the  conditions  of  Fig. 

1,  we  have  C(R)  ■  0.038  -  cst  and  eq.  5  can  be  used  with  no  further 
approximation  to  calculate  ag(R).  The  1-mrad  f ield-of-view  (half 
angle)  signal  S(R)  is  used  for  this  application.  This  field  of  view  is 
sufficiently  small  to  minimize  the  multiscattering  contributions  and  it 
simulates  well  practical  applications. 

Figure  6  compares  the  true  profile  of  the  extinction  coefficient 
with  the  solutions  obtained  by  specifying  the  value  of  af  at  the  near¬ 
end  boundary  of  the  cloud.  Three  solutions  were  calculated  for  * 
0.9,  1.0  and  1.1  km-*  respectively,  the  true  value  being  1.0  km-1.  As 
shown,  the  solutions  for  *  1.0  and  1.1  km-1  diverge  to  infinity 
while  that  for  ■  0.9  km-1  quickly  falls  off  to  zero.  Even  the  solu¬ 
tion  for  the  true  value  of  1.0  km-1  diverges  because  the  1-mrad 
simulated  lidar  signal  contains  enough  multiscattering  contributions  to 
cause  irrecoverable  instability.  Therefore,  Fig.  6  shows  that  the 
solutions  based  on  the  near-end  boundary  condition  are  totally  defi¬ 
cient  and  unacceptable  for  conditions  represented  by  the  cloud  defined 
in  Fig.  1.  By  comparison.  Fig.  7  shows  the  solutions  for  the  same 
values  of  a^but  specified  at  the  far-end  boundary  of  the  cloud  where 
the  true  ag  is  also  1.0  km-1.  All  three  solutions  are  acceptable;  the 
original  errors  in  are  damped  and  the  three  curves  converge  to  a 
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single  one  as  the  front  end  of  the  cloud  Is  approached.  The  small 
differences  with  the  true  profile  come  from  the  multiscattering  contri¬ 
butions  which  are  not  taken  into  account  by  the  single-scattering  lidar 
equation.  In  summary,  the  results  illustrated  in  Figs.  6  and  7  only 
confirm  previous  findings  (Refs.  1  and  5)  that  the  stability  of  the 
single-scattering  lidar  equation  requires  that  the  boundary  condition 
be  given  or  specified  at  the  far  end  of  the  lidar  range  where  the 
signal  S(R)  is  generally  smallest. 

Figure  8  shows  the  solutions  obtained  for  the  same  cloud  and 
with  the  same  method  as  in  Fig.  7  but  given  large  errors  in  a^,  namely 
*  0.1  and  3.0  km-*  compared  with  the  true  value  of  1.0  km-1.  The 
solutions  eventually  converge  as  expected,  but  the  errors  remain  large 
over  a  nonnegligible  portion  of  the  cloud.  This  is  especially  true  in 
the  case  of  a ^  **  0.1  km-1  where  the  differences  normalized  to  the  cloud 
average  extinction  coefficient  of  3.5  km-*  can  grew  as  high  as  80% 
before  being  damped.  Therefore,  a  reasonable  estimate  of  on  the  far 
side  of  the  cloud  Is  an  important  requirement  even  though  large  errors 
do  not  cause  irrecoverable  instabilities.  Good  accuracy  In  estimating 
or  independently  measuring  is  difficult  to  achieve  since  the  far- 
side  boundary  Is  shielded  from  the  lidar  station  by  the  cloud  itself. 
For  clouds  of  arbitrary  extinction  profiles,  there  Is  yet  no  systematic 
method  of  deriving  from  lidar  data. 

The  example  of  Fig.  8  as  well  as  many  other  calculations  show 
that  the  solutions  converge  more  rapidly  If  is  overestimated  rather 
than  underestimated.  For  instance,  the  solution  based  on  *  3.0  km-* 
becomes  acceptable  at  about  the  2/3  point  of  the  cloud  (measured  from 
the  front  end),  while  this  occurs  only  at  the  1/3  mark  for  the  solution 
based  on  af  ■  0.1  km-1.  This  is  in  agreement  with  the  theoretical 
analysis  of  Ref.  5  (eq.  15)  which  shows  that  the  relative  error  on  the 
calculated  extinction  coefficient  is  inversely  proportional  to  the 
boundary  value. 
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As  discussed  earlier,  a  fundamental  characteristic  of  the 
single-scattering  lidar  inversion  problem  described  by  eqs.  1-5  is  that 
a  relation  between  the  extinction  and  backscattering  coefficients  must 
be  assumed*  We  have  seen  that  the  simplification  C(R)  -  cst  can  be 
made  if  the  composition,  geometry  and  size  distribution  of  the  aerosol 
particles  remain  invariant.  We  will  now  investigate  if  this  simplifi¬ 
cation  can  be  applied  to  more  realistic  clouds  of  variable  C(R).  The 
cloud  defined  in  Fig.  2  has  the  same  extinction  profile  as  In  Fig.  1 
but  it  differs  by  a  smoothly  varying  function  C(R).  The  solutions 
obtained  for  this  cloud  under  the  assumption  of  a  constant  C(R)  are 
plotted  in  Fig.  9.  The  boundary  values  are  as  in  Fig.  8,  i.e.  ■ 
0.1,  1.0  and  3.0  km-1  specified  at  the  far-end  boundary  of  the  cloud. 
The  solutions  are  not  very  much  different  from  those  of  Fig.  8  for  a 
truly  constant  C(R).  They  tend  to  be  greater  in  regions  where  C(R)  is 
greater  than  its  value  at  the  boundary  point  and  vice  versa,  which  is 
in  agreement  with  the  analytic  expressions  (eqs.  12  and  15)  derived  in 
Ref.  5.  The  symmetry  of  the  calculated  extinction  distribution  is  thus 
skewed,  even  for  the  exact  boundary  value.  However,  the  largest  errors 
are  again  associated  with  the  boundary-value  errors. 

The  solutions  for  the  more  complex  case  of  Fig.  3,  where  regions 

of  constant  but  different  C(R)  ■  8/a  are  juxtaposed,  are  shown  in  Fig. 

e 

10.  For  the  exact  a^,  the  solution  remains  exact  so  long  as  C(R)  does 
not  change  from  its  value  at  the  far-end  boundary.  With  the  step 
changes  in  C(R),  the  computed  extinction  coefficient  becomes  In  error 
by  an  amount  that  appears  to  be  proportional  to  the  magnitude  of  the 
step.  If  C(R)  increases,  the  extinction  coefficient  becomes  greater 
than  the  true  value  and  vice  versa.  For  the  conditions  of  Fig.  10,  the 
differences  resulting  from  the  assumption  of  constant  C(R)  can  become 
as  high  as  2/3  times  the  local  true  extinction  coefficient.  It  is 
emphasized,  however,  that  nowhere  in  Figs.  9  and  10,  nor  in  other  cal¬ 
culations,  did  the  assumption  of  constant  C(R)  lead  to  catastrophic 
results  as  in  Fig.  6  for  the  near-end  solution. 
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The  solutions  for  the  clouds  defined  in  Figs.  2  and  3  were 
recalculated  after  substitution  in  eq.  5  of  their  respective  true  func¬ 
tion  C(R).  The  results  are  plotted  in  Figs.  11  and  12  respectively. 

The  solutions  for  the  true  boundary  value  reproduce  almost  exactly 
the  true  extinction  profiles.  The  solutions  for  the  erroneous 
remain  in  error  in  the  far-end  region  of  the  cloud  but  gradually  tend 
toward  the  true  profile  with  decreasing  range.  As  previously  noted, 
the  errors  are  less  if  the  boundary  value  is  overestimated.  Finally, 
it  can  be  observed  that  the  solutions  of  Fig.  11  are  very  close  to 
those  of  Fig.  8;  the  true  extinction  profiles  are  indeed  the  same  in 
both  cases,  only  the  functions  C(R)  are  different.  This  result  con¬ 
firms  that  the  single-scattering  lidar  solution  gives  the  correct 
extinction  coefficient  if  C(R)  is  known. 

The  solutions  of  Figs.  11  and  12  show  that  the  knowledge  of  C(R) 
allows  the  determination  of  the  extinction  coefficient  profile  with  a 
precision  dictated  only  by  that  of  its  boundary  value  at  the  far  side 
of  the  cloud  and,  of  course,  by  the  experimental  errors.  However,  they 
give  no  clue  as  to  how  this  Input  can  be  obtained  In  practice.  The 
lidar  signal  alone  Is  not  sufficient.  Additional  information  on  the 
nature  of  the  aerosol  particles  is  required,  for  example  In  the  form 
of  an  empirical  correlation  between  p  and  ag  as  suggested  by  Klett 
(Ref.  2). 

The  solutions  of  Figs.  8  and  11  calculated  with  the  exact 
boundary  value  and  the  prescribed  C(R)  functions  differ  slightly 
although  the  true  extinction  profile  is  the  same  in  both  cases.  The 
reason  Is  multiscattering  contributions.  These  contributions  are  more 
important  in  the  case  of  Fig.  8  because  the  forward-scattering  peak  is 
narrower  for  the  99%  relative  humidity  maritime  aerosol  (Fig.  8)  than 
for  the  mixture  of  rural,  urban  and  maritime  aerosols  at  702  relative 
humidity  (Fig.  11).  Indeed,  a  narrower  forward  peak  means  that  more 
forward-scattered  radiation  remains  in  the  receiver  field  of  view.  As 
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illustrated  in  Figs.  8  and  11,  the  multiscattering  contributions  have 
the  effect  of  reducing  the  extinction  coefficient  calculated  with  the 
single-scattering  lidar  equation.  This  is  easily  understood  from  eq. 
1.  The  increase  of  the  lidar  signal  by  multiscattering  returns  is 
interpreted  by  eq.  1  as  a  drop  in  the  extinction  coefficient  since  the 
exponential  is  a  faster  function  than  the  linear  backscatter  term 
(proportional  to  a  according  to  eq.  4). 


The  multiscattering  effect  depends  on  several  factors.  First, 
there  is  the  nature  of  the  aerosols  which  explains,  as  we  have  just 
seen,  the  differences  between  Figs.  8  and  11.  The  important  property 
is  the  width  of  the  forward  scattering  peak;  the  narrower  the  width, 
the  greater  the  probability  of  recording  multiscattering  contributions 
by  a  receiver  of  given  field  of  view.  Second,  there  is  the  cloud  opti¬ 
cal  depth.  For  a  given  aerosol,  the  multiscattering  returns  increase 
with  the  optical  depth.  For  example.  Fig.  13  shows  what  happens  if  the 
aerosol  concentration  of  the  cloud  defined  in  Fig.  1  is  doubled.  Com¬ 
paring  the  results  depicted  in  Figs.  8  and  13,  we  observe  that  the 
single-scattering  solution  underestimates  substantially  the  extinction 
coefficient  at  greater  optical  depths.  Finally,  an  obvious  factor  is 
the  receiver  field  of  view.  Figure  14  contrasts  the  solutions  derived 
from  lidar  signals  obtained  from  the  same  cloud  with  receiver  fields  of 
view  of  1,  5  and  20  mrad.  The  predicted  extinction  coefficient 
decreases  with  widening  field  of  view,  which  is  in  agreement  with  the 
preceding  observations  on  the  effect  of  increasing  multiscattering 
contributions.  In  practice,  different  cloud-to-receiver  distances  will 
produce  the  same  result-  This  is  illustrated  in  Fig.  15  where  the 
solutions  for  the  1-mrad  lidar  signals  are  presented  for  cloud-to- 
receiver  distances  of  0.2,  1.0  and  3.0  km. 


Based  on  the  preceding  results,  the  following  observations  can 
be  made  on  the  use  of  the  single-scattering  lidar  equation  to  determine 
the  aerosol  extinction  coefficient: 
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a)  At  extinction  coefficients  of  the  order  of  1  km""1  or 
greater,  the  near-end  boundary  value  inversion  algorithm  is 
unstable  to  small  errors  in  the  assumed  boundary  value  and 
to  small,  but  unavoidable,  multiscattering  contributions. 

The  same  is  probably  true  for  experimental  errors.  The 
instabilities  are  irrecoverable  and  make  the  method  useless 
in  most  applications  involving  clouds. 

b)  The  far-end  boundary  value  inversion  algorithm  is  stable  to 
any  perturbations  resulting  from  boundary  value  errors, 
uncertainties  in  the  backscatter-to-extinction  ratio,  and/or 
multiscattering  contributions.  The  remaining  comments  con¬ 
cern  solutions  obtained  with  this  method. 

c)  It  is  the  boundary  value  errors  that  most  affect  the  calcu¬ 
lated  extinction  coefficient.  Although  the  solutions  even¬ 
tually  converge  to  the  true  profiles,  the  convergence  rate 
depends  on  the  initial  gap,  and  the  residual  errors  may  be 
large.  For  all  cases  studied,  the  convergence  rate  is 
faster  if  the  assumed  far-end  extinction  coefficient  Is 
overestimated.  The  determination  of  the  boundary  value  is  a 
serious  problem,  and  the  information  contained  in  the 
single-scattering  lidar  signal  is  not  sufficient  to  permit  a 
reliable  estimate  for  arbitrary  extinction  profiles. 

d)  The  errors  resulting  from  the  assumption  of  a  constant 
profile  of  the  backscatter-to-extinction  ratio  (3/ ag  are, 
surprisingly,  within  acceptable  limits.  For  step  changes  in 
p/ag  of  the  order  of  4,  errors  in  the  local  extinction 
coefficient  may  exceed  100%,  but  the  general  shape  of  the 
profile  is  obtained  and  the  average  aerosol  extinction 
coefficient  over  the  cloud  depth  is  within  10-20%  of  the 
true  value.  If  the  function  0/afi  is  known,  the  solution 
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converges  to  the  true  extinction  profile.  For  practical 

applications,  additional  and  independent  information  on  the 

nature  of  the  aerosol  particles  is  required  to  determine  or 

estimate  8/a  . 

r  e 

e)  If  the  receiver  field  of  view  is  chosen  to  stare  only  at  the 
volume  probed  by  the  unscattered  transmitter  beam  and  if 
that  beam  is  reasonably  collimated,  the  solution  errors 
caused  by  multiscattering  contributions  remain  below  a  few 
percent  up  to  an  optical  depth  between  3  and  A.  The  degra¬ 
dation  beyond  this  limit  is  quite  rapid,  but  it  is  question¬ 
able  whether  lidar  measurements  remain  reliable  for  optical 
depths  greater  than  3-4. 
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FIGURE  6 


FIGURE  7 


Range  (km) 

-  Single-scattering  lidar  solutions  for  the  cloud  defined  in 
Fig.  1  calculated  with  the  near-end  boundary  value  algo¬ 
rithm.  The  solutions  were  obtained  with  extinction  coeffi¬ 
cient  boundary  values  (B.V.)  of  0.9,  1.0  and  1.1  km-1  for  a 
true  value  of  1.0  km-1. 


RANGE  (km) 

-  Single-scattering  lidar  solutions  for  the  cloud  defined  in 
Fig.  1  calculated  with  the  far-end  boundary  value  algorithm 
The  solutions  were  obtained  with  extinction  coefficient 
boundary  values  (B.V.)  of  0.9,  1.0  and  1.1  km-1  for  a  true 
value  of  1.0  km-1. 


UNCLASSIFIED 

20 


RANGE  (km) 


FIGURE  10 


Same  as  Fig. 


9  but  for  the  cloud  defined  In  Fig.  3. 


RANGE  (km) 


FIGURE  11  -  Same  as  Fig.  9  except  that  the  solutions  were  calculated 
with  the  true  p/ae  profile. 
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RANGE  (km) 


FIGURE  12  -  Same  as  Fig.  10  except  that  the  solutions  were  calculated 
with  the  true  {3/ae  profile. 


RANGE  (km) 


FIGURE  13  -  Single-scattering  lidar  solutions  calculated  with  the 

far-end  boundary  value  algorithm  for  the  cloud  defined  in 
Fig.  1  except  that  the  extinction  profile  was  doubled  to 
exemplify  the  multiscattering  effects.  The  solutions  were 
obtained  with  boundary  values  (B.V.)  of  0.2,  2.0  and  6.0 
km-^  for  a  true  value  of  2.0  km-1. 
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RANGE  (km) 

FIGURE  14  -  Single-scattering  lidar  solutions  calculated  with  the 

far-end  boundary  value  algorithm  for  the  lidar  signals  from 
the  cloud  defined  in  Fig.  1  generated  with  fields  of  view 
(FOV)  of  1,  5  and  20  mrad.  The  true  boundary  value  of  1.0 
km-1  was  used  for  all  three  solutions. 


RANGE  IN  CLOUD  (km) 

FIGURE  15  -  Single-scattering  lidar  solutions  calculated  with  the 

far-end  boundary  value  algorithm  for  the  lidar  signals  from 
the  cloud  defined  in  Fig.  1  generated  with  a  field  of  view 
of  1  mrad  but  cloud-to-receiver  distances  (cloud-rec' r)  of 
0.2,  1.0  and  3.0  km.  The  true  boundary  value  of  1.0  km-i 
was  used  for  al  solutions. 


UNCLASSIFIED 

23 


A.O  MULTI SCATTERING  INVERSION  METHOD 


A.l  Solution  Method 


The  proposed  lidar  solution  method  is  derived  from  the  multi¬ 
scattering  propagation  model  described  in  Ref.  10.  The  aim  is  to  make 
use  of  the  information  carried  by  the  higher  order  scattering  events 
occurring  on  both  transmission  legs  to  and  from  the  probed  volume. 

From  the  solution  of  Ref.  10  for  the  backscattered  radiant  intensity, 
the  power  detected  by  a  lidar  receiver  that  shares  the  same  axis  as  the 
transmitter  is  given  by 

PA  R 

P(R,u>)  *  — 2—  a”(R)  exp [-  2  /(a^  +  afl  +  a“)dr] 

xw2  o 

o 

x  |  exp[-  J  a+  dr]  +  V(R,(j)  /  — — -  a+(r*)exp[-  J  orkdr]  j  [6] 

l  K(R)  o  S  o  M(R.r’)  S  o  8  \ 

where  Pq  Is  the  transmitter  power,  wq  is  the  zero-range  radius  of  the 

transmitter  beam,  A  Is  the  collecting  area  of  the  receiver,  R  Is  the 

range,  w  is  the  receiver  field  of  view,  and  a  is  the  molecular  absorp- 

m 

tion  coefficient.  The  aerosol  scattering  and  absorption  properties  are 
given  by  the  backscattering  coefficient 

x  2x  ®  dQ 

a"  *  /  d6  sin0  /  d(J>  /  dp  ■  j  P  xp2  (9,<t>,p),  [7] 

x/2  o  o  p 

the  forward  scattering  coefficient 

x/2  2x  ®  .  dQ 

atg  "  /  d0  slnQ  /  dip  J  dp  P  xp2  ( 0 , <J>,  p) ,  [8] 

o  o  o  ^ 


and  the  aerosol  absorption  coefficient 
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The  remaining  functions  U(R,u)  and  V(R,w)  are  f ield-of-view 
reduction  factors  given  in  Ref.  10.  They  depend  on  oj  and  the  scatter¬ 
ing  properties  of  the  aerosols.  They  are  complicated  functions  of 
these  arguments  but  only  their  asymptotic  values  will  be  used  here. 
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For  optical  and  infrared  lidar  applications  in  the  atmosphere, 
the  scattering  phase  function  p( 9)  is  much  more  sharply  peaked  In  the 
forward  than  in  the  backward  direction.  It  follows  from  eqs.  12  and  13 
that  D+/(r-r')  <<  D"/(R-r')  and  we  can  neglect  the  third  term  on  the 
right-hand  side  of  eq.  11.  Hence, 

M(R.r’)  =  K(R),  [14] 

and  eq.  6  becomes 

PA  R 

P(R,cj)  - - - - a~(R)  exp[-  2  /  (a  +  a  +  a~)dr] 

o  «•  /  _  \  s  L  m  a  s  J 

nw^K(R)  o 

o 


Equation  15  is  the  proposed  modified  lidar  equation  that  takes 
into  account  multiple  scatterings.  It  is  instructive  to  consider  two 
asymptotic  forms  of  eq.  15.  First,  we  examine  the  case  of  an  open 
receiver,  I.e.  a  field  of  view  (half  angle)  w  equal  to  it/2 .  By  defini¬ 
tion  (Ref.  10),  U  and  V  tend  to  unity  in  that  limit  and  eq.  15  becomes 

PA  R 

P(R,x/2)  *  - - - a~(R)  exp[-  2  /  (a  +  a  +  a~)dr].  [16] 

*w2K(R)  5  FL  J  ■  a  s  J  L  j 

Second,  we  consider  the  situation  where  the  receiver  field  of  view 
8 tares  at  the  volume  probed  by  the  unscattered  transmitter  beam  only. 

In  other  words,  the  field  of  view  is  chosen  equal  or  just  slightly 
greater  than  the  transmitter  beam  divergence  assumed  to  be  small.  For 
that  condition,  it  can  be  shown  from  the  expressions  of  Ref.  10  for  U 
and  V  that 

lim  V(R, to)  =  0,  [17] 

0)  -*■  w  L  J 


[15] 


IV  A. 

U(R,w)  exp[-  /  a+dr]  +  V(R,oj)[1  -  exp(-  /  a+dr)  ] 


Mi 


JIT, 


Si 


»  *»  -  v*  ’ 
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lim  U(R,w)  =  exp[-  /  a+dr],  [l8] 

where  is  equal  to  the  transmitter  divergence.  Physically,  eqs.  17 

and  18  state  that  mjltiscattering  contributions  are  not  significant  at 

such  small  fields  of  view:  the  fraction  V  of  the  multiscatter  signal 

that  would  originate  from  the  forward-scattered  beam  is  approximately 

zero  since  the  receiver  field  of  view  sees  only  a  small  region  of  that 

broadened  beam,  and  the  fraction  U  arising  from  the  unscattered  beam 

R 

falls  off  as  exp[-  /  a+dr]  since  the  radiation  scattered  on  the  trans- 

o 

mission  path  back  to  the  receiver  is  collected  in  a  negligible  amount 

because  of  the  small  field  of  view.  The  approximations  of  eqs.  17  and 

R 

18  begin  to  fail  at  optical  scattering  depths  J  a+dr  beyond  3-4. 

o 

Substituting  eqs.  17  and  18  for  U  and  V  in  eq.  15,  we  obtain 

PA  R 

P(R.«t)  "  — r~ -  “g(R>  e*P[“  2  /  (am  +  “a  +  “s  +  °£)drl*  f193 

mw2K(R)  o 

o 

Except  for  slight  differences  in  the  formulations  of  the  backscatter 

coefficient  a-  and  geometric  factor  K(R),  eq.  19  is  similar  to  the 
s 

single-scattering  lidar  equation.  Both  give  approximately  the  same 
answers  for  P  as  demonstrated  by  the  results  of  Ref.  10. 


If  we  ratio  eqs.  16  and  19,  we  obtain 


P(R,w/2)  _ 

iw,wtr 


R 

exp [  2  /  a+d z  ] . 
s 
o 


[20] 


Equation  20  is  at  the  origin  of  the  proposed  multiscattering  lidar 
inversion  method.  If  we  could  record  simultaneously  the  lidar  returns 
at  both  fields  of  view  u> and  n/2,  we  could  easily  determine  the 
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forward-scattering  coefficient  from  eq.  20.  The  technique  would 
require  no  instrument  calibration,  no  relationship  between  backscatter 
and  extinction,  and  no  boundary  value.  It  could  thus  be  free  of  the 
most  serious  difficulties  associated  with  the  solution  of  the  single¬ 
scattering  lidar  equation.  The  parameter  determined  would  be  the 
forward-scattering  coefficient  which  is,  within  a  few  percent,  equal  to 
the  albedo  times  the  extinction  coefficient. 

The  technique  suggested  by  eq.  20  is  straightforward,  but  lidar 
measurements  at  wide  fields  of  view,  as  required  for  P(R,it/2),  are  not 
practical  because  of  signal-to-noise  problems.  Nonetheless,  the 
ratioing  principle  can  still  be  applied  to  simultaneous  measurements  at 
smaller  fields  of  view  and  still  carry  the  advantages  of  no  require¬ 
ments  for  a  boundary  value  and  for  a  backscatter-to-extinction  rela¬ 
tion.  Equation  20  would  have  to  be  replaced  by  an  expression  involving 
the  functions  U  and  V.  Since  U  and  V  are  complicated  functions  of 
their  arguments,  they  do  not  provide,  in  general,  a  mathematical  form 
useful  for  inversion  algorithms.  However,  calculation  results  and  the 
asymptotic  form  of  eq.  20  suggest  the  following  simple  empirical 
formula: 

iSr  ln  "  2  f(u,U)t}  °s(R)’  t21 1 

where  u  is  a  receiver  field  of  view  greater  than  u)t,  and  f(u),u>t)  is  a 
function  of  w  and  to  be  determined  empirically.  f(o>,ut)  is  not  a 
universal  function,  it  depends  on  the  aerosol  properties. 

By  definition,  Is  small  and  approximately  equal  to  the  trans¬ 
mitter  beam  divergence.  On  the  other  hand,  the  field  of  view  to  must  be 
chosen  large  enough  to  make  f(w,ut)  measurable  but  as  small  as  possible 
to  minimize  the  noise  problems.  Calculation  results  obtained  under 
these  conditions  suggest  the  following  empirical  form: 
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f(u),u>t)  -  (x/b)  . 


/V  +  f'tJk 
.  o 

X  “  In  - x -  . 

/TT”  Ao.  R 
o  t 


for  u  <  u  «  n/2,  where  Wq  is  the  radius  of  the  receiver  aperture. 

The  parameter  b  appears  to  be  a  system  constant,  but  k  depends  on  the 
aerosol  properties  and  can  vary  with  range  R.  Figure  16  shows  hew  eq. 
22  compares  with  calculated  lidar  signals  for  two  types  of  aerosols. 

The  power  law  function  provides  a  very  good  fit  up  to  fields  of  view 
(half  angle)  of  the  order  of  100  mrad.  The  aerosols  chosen  represent 
conditions  that  are  as  far  apart  as  possible  for  the  AFGL  generic 
models  at  the  1.06-pm  wavelength.  The  backscatter-to-extinction  ratio 
and  albedo  are  0.038  and  0.998,  and  0.010  and  0.634  for  the  maritime 
and  urban  aerosols  used  for  the  calculations  illustrated  in  Fig.  16, 
respectively.  Calculations  with  different  clouds  produced  similar  fits 
with  local  discrepancies  less  than  5%  for  fields  of  view  (half  angle) 
smaller  than  about  100  mrad.  The  empirical  model  defined  by  eqs.  22-23 
constitutes  one  convenient  choice  but  some  other  functional  dependence 
could  probably  give  similar  or  even  better  results. 


Equations  21-23  form  the  basis  of  a  very  simple  lidar  inversion 
technique.  If  the  lidar  return  is  recorded  simultaneously  at  three 


fields  of  view  wj  and  u)2»  the  exponent  k  can  be  calculated  with  the 


formula 


lnjjffi-  ln[p(R,uI)/P(R,u)t)^r  ln[P(R,h)2)/P(R,h>t)] 

ln[xi/x2] 


.  [2*] 


where  Xj  and  X2  are  obtained  by  substitution  of  and  u>2  for  co  in  eq. 
23.  In  practice,  the  ratios  P(R,oJ1)/P(R,wt )  and  P(R,  co2)/P(R,wt) 


m 


vms&mm,' 


/.  •<  *  ■ 


■C-V 
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should  be  smoothed  to  minimize  the  fluctuations  resulting  from  the 
finite  difference  approximations  to  the  derivatives  of  eq.  24. 

Once  k(R)  is  determined,  the  forward-scattering  coefficient  a+ 

s 

is  calculated  with  the  following  formula  derived  from  eq.  21: 

<z+(R)  -  I  (b/X)k  ^  ln[P(R,w)/P(Riu>t)],  [25] 

for  w  *  uij  and/or  io2t  and  x  is  given  by  eq.  23.  The  parameter  b  is  a 
system  constant  and  must  be  determined  by  calibration  or  other  means. 
The  results  plotted  in  Fig.  16  show  that  different  aerosols  give  dif¬ 
ferent  exponents  k.  Since  k  is  computed  at  each  range  R,  the  method  is 
therefore  applicable  to  clouds  of  varying  composition. 

Equations  24  and  25  define  the  proposed  multiscattering  lidar 
inversion  algorithm.  It  is  based  on  the  simultaneous  measurement  of 
the  lidar  returns  at  three  fields  of  view.  One  field  of  view  must  be 
chosen  equal  to  or  just  slightly  above  the  transmitter  divergence.  The 
other  two  are  greater  but  Just  enough  to  give  measurable  differences. 
Since  only  ratios  are  used  in  the  algorithm,  no  individual  signal  cali¬ 
bration  is  necessary.  The  ratios  are  simply  adjusted  to  tend  to  unity 
at  the  front  end  of  the  cloud  or  zero  optical  depth.  The  solution  for 
the  forward  scattering  coefficient  a*  is  calculated  by  straightforward 
substitution  in  the  algebraic  formulas  23-25.  No  backscatter-to- 
extinction  relation  and  no  boundary  value  are  required.  Equation  25 
shows  that  the  proposed  inversion  method  is  very  similar  to  the  simple 
slope  method  but  applied  to  the  lidar  ratio  P(R,u))/P(R,u>t )  • 
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4.2  Results 


The  proposed  multiscattering  lidar  inversion  technique  is 
applied  to  the  signals  generated  for  the  clouds  defined  in  Figs.  1-4. 
Examples  of  lidar  ratios  used  for  the  calculations  are  shown  in  Fig. 

17,  they  were  obtained  for  the  cloud  profiles  of  Fig.  3.  These  results 
illustrate  the  magnitude  of  the  signal  differences  that  can  be  expected 
in  atmospheric-like  aerosols  for  fields  of  view  (half  angle)  of  1,  20, 
50  and  100  mrad.  Such  differences  are  measurable  with  acceptable 
signal-to-noise  ratios.  Near  simultaneous  measurements  of  lidar 
returns  at  1  and  3  mrad  off  marine  stratus  clouds  with  a  480-m  altitude 
cloud  base  (Ref. 6)  showed  a  signal  ratio  that  can  reach  a  value  of  2. 
Therefore,  the  calculated  ratios  of  order  10  or  less  plotted  in  Fig.  17 
appear  reasonable  considering  the  greater  differences  in  fields  of 
view. 


The  results  for  the  cloud  of  uniform  composition  defined  In 

Fig.  1  are  plotted  in  Fig.  18.  They  were  calculated  with  and  u>2 

equal  to  1,  20  and  50  mrad.  As  shown,  the  solution  fits  the  true  a+ 

profile  very  well,  absolute  local  errors  are  less  than  0.1  km-1  or 

about  3 %  of  the  average  true  a+.  The  cloud-integrated  a+,  or  forward- 

s  s 

scattering  optical  depth,  agrees  with  the  true  value  within  1%. 

The  results  for  the  cloud  defined  in  Fig.  2  are  shown  in  Fig. 

19.  Despite  important  variations  in  albedo  and  backscatter-to- 
extinction  ratio,  the  agreement  with  the  true  a+  is  still  acceptable, 
local  errors  are  less  than  10%.  The  solutions  for  the  more  complex 
cloud  structures  illustrated  in  Figs.  3  and  4  are  given  in  Figs.  20  and 
21,  respectively.  Again,  the  agreement  is  very  good  with  some  excep¬ 
tions  at  points  where  the  albedo  and  p/a  change  abruptly.  The  curaula- 
tive  effects  of  the  solution  errors  on  the  calculated  cloud  optical 
depth  are  less  than  10-15%.  In  all  cases,  the  true  profiles  are  well 
reproduced  by  the  solutions. 
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As  discussed  previously,  the  proposed  nultiscattering  technique 
requires  that  the  lidar  return  at  the  smallest  field  of  view  u>t  be 
dominated  by  single  backscatterings .  However,  at  very  large  optical 
depths,  the  lidar  signal  at  u>  may  be  contaminated  by  multiple  scatter¬ 
ings.  The  solution  obtained  for  the  cloud  defined  in  Fig.  1  with  the 
extinction  profile  multiplied  by  a  factor  of  2  is  drawn  in  Fig.  22. 

The  calculations  were  performed  with  u>t ,  and  oj2  equal  to  1,  20  and 
50  mrad  as  for  the  examples  given  in  Figs.  18-21.  It  is  shown  that  the 
solution  falls  to  reach  the  high  scattering  coefficient  value  at  the 
cloud  center  and  that  it  remains  below  the  true  a*  profile  from  that 
point  on.  However,  the  solution  does  not  diverge  and  the  profile  has 
approximately  the  right  shape. 

The  results  of  Fig.  22  Indicate  that  the  multiscattering  lidar 
solution  method  becomes  subject  to  important  errors  for  optical  depths 
greater  than  about  4.  Although  this  limit  is  expected  to  vary  with  the 
nature  of  the  aerosols  and  the  profile  of  the  extinction  coefficient, 
it  is  not  expected  that  the  method  defined  by  eqs.  23-25  will  remain 
accurate  beyond  optical  depths  of  the  order  of  4.  This  is  similar  to 
the  limit  found  In  the  preceding  chapter  for  the  single-scattering 
inversion  method  under  optimal  conditions  of  known  boundary  value  and 
p/cte  profile.  This  is  not  surprising  since  the  requirement  on  P(R,wt>, 
namely  that  it  Is  not  significantly  affected  by  multiple  scatterings, 
is  the  same  as  for  the  validity  of  the  single-scattering  lidar  equation, 
A  more  sophisticated  algorithm  involving  measurements  at  more  than 
three  fields  of  view  could  possibly  be  defined  to  extend  the  range  of 
applicability.  However,  it  is  questionable  whether  meaningful  lidar 
measurements  can  be  carried  out  at  optical  depths  greater  than  4. 
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FIGURE  16  -  Function  f(<j,  ii>t)  defined  In  eq.  21  calculated  for  the 

generic  models  of  maritime  aerosol  at  99 2  relative  humidity 
and  urban  aerosol  at  702  relative  humidity  (Ref.  11).  The 
wavelength  Is  1.06  pm.  The  symbols  are  calculations  with 
the  propagation  model  of  Ref.  10  and  the  continuous  curve 
represents  empirical  fit  of  eq.  22.  k  ■  1.34  for  maritime 
aerosol  and  2.42  for  urban  aerosol;  b  »  2.6  in  both  cases. 


FIGURE  17  -  Examples  of  ratios  of  Ildar  signals  at  different  fields  of 
view  calculated  from  the  cloud  defined  In  Fig.  3. 
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5.0  SUMMARY  AND  CONCLUSIONS 

This  report  has  described  a  new  lidar  inversion  method  for  the 
determination  of  the  aerosol  scattering  coefficient.  Its  principal 
characteristic  is  the  use  of  the  information  contained  in  the  multi¬ 
scattering  contributions  of  the  measured  lidar  returns.  The  method 
requires  neither  a  boundary  value  nor  a  backscatter-to-extinction  rela¬ 
tion.  Performance  proved  to  be  more  consistent  than  with  the  standard 
single-scattering  lidar  method  for  the  same  simulated  applications. 

The  simulated  lidar  signals  were  generated  with  the  multiscat¬ 
tering  propagation  model  of  Ref.  9.  Inhomogeneous  clouds  were  con¬ 
structed  from  the  generic  maritime,  rural  and  urban  aerosol  models 
developed  by  the  Air  Force  Geophysics  Laboratory  (Ref.  11).  These 
models  were  designed  from  measurements  and  they  represent  adequately  a 
wide  variety  of  atmospheric  conditions.  Range-resolved  lidar  signals 
were  calculated  as  functions  of  the  receiver  field  of  view. 

The  slngle-6cattering  lidar  Inversion  algorithm  based  on  the 
specification  of  the  boundary  value  at  the  front  end  of  the  cloud  was 
found  to  be  inadequate  under  all  conditions  studied.  The  solutions  for 
the  aerosol  extinction  coefficient  either  diverged  to  infinity  or 
quickly  fell  off  to  zero.  On  the  other  hand,  the  inversion  method  that 
integrates  from  the  backside  of  the  cloud  onward  remained  6table  In  all 
cases.  However,  the  results  showed  that  the  calculated  extinction 
coefficient  can  be  seriously  affected  by  uncertainties  in  the  boundary 
value  and  in  the  backscatter-to-extinction  ratio  for  clouds  of  inhomo¬ 
geneous  composition. 

The  proposed  multiscattering  lidar  inversion  method  requires 
that  the  range-gated  lidar  signals  be  measured  simultaneously  at  three 
fields  of  view.  One  field  of  view  is  small  and  corresponds  to  the 
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standard  design  where  the  receiver  stares  only  at  the  volume  probed  by 
the  unscattered  transmitter  beam;  the  other  two  must  be  chosen  just 
wide  enough  for  the  differences  due  to  multiple  scatterings  to  be  mea¬ 
surable  with  sufficient  accuracy.  The  aerosol  scattering  phase  func¬ 
tion  must  be  sufficiently  peaked  forward  for  these  fields  of  view  to  be 
kept  within  practical  limits.  The  inputs  to  the  solution  method  are 
the  ratios  of  the  lidar  signals  at  the  different  fields  of  view.  The 
algorithm  Involves  only  algebraic  equations  and  is  very  similar  to  the 
well-known  slope  method.  It  requires  the  calibration  of  one  system 
constant:  the  parameter  b  in  eq.  22.  The  calculation  results  fit  very 

well  the  true  profiles  in  all  cases.  The  agreement  was  consistently 
better  than  for  the  single-scattering  solutions  except  where  the  exact 
boundary  value  and  the  exact  backscatter-to-extinction  profiles  were 
given,  in  which  case  performances  are  similar.  The  algorithm  could 
still  be  Improved,  particularly  the  determination  of  the  f ield-of-view 
function  f(u),ut)  defined  In  eq.  21. 

The  method  solves  for  the  aerosol  forward-scattering  coefficient 
defined  by  eq.  7  and  not  for  the  extinction  coefficient,  but  the  two 
are  related  through  a  proportionality  constant.  For  atmospheric 
aerosols  and  operating  wavelengths  in  the  visible  and  Infrared,  the 
proportionality  constant  is  approximately  equal  to  the  aerosol  single- 
scattering  albedo;  It  is  actually  smaller,  by  5-10%  at  the  most. 

Except  for  highly  absorbing  aerosols,  the  differences  between  the 
forward-scattering  and  extinction  coefficients  are  probably  within  the 
experimental  errors. 

With  the  standard  single-scattering  solution  method,  the  tacti¬ 
cal  use  of  lidar  for  the  determination  of  the  atmospheric  aerosol 
optical  and  Infrared  parameters  is  seriously  hindered  by  the  need  for 
advanced  knowledge  of  a  boundary  value  at  the  far  end  of  the  lidar 
range  and  of  the  profile  of  the  backscatter-to-extinction  ratio.  The 
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proposed  multiscattering  lidar  inversion  technique  eliminates  these 
difficulties.  The  additional  information  is  extracted  from  the  multi- 
scattering  contributions  through  simultaneous  measurements  of  the  lidar 
returns  at  three  fields  of  view.  The  aerosol  forward-scattering  coef¬ 
ficient  is  thus  determined  from  lidar  data  only,  independently  of  a 
boundary  value  and  of  the  backscatter-to-extinction  profile.  Under 
such  conditions,  the  use  of  lidar  as  a  tactical  aid  to  provide  reliable 
aerosol  data  would  be  feasible.  Of  course,  the  technique  needs  experi¬ 
mental  validation.  Calculations  have  shewn  that  receiver  fields  of 
view  in  the  range  of  1,  20  and  50  mrad  (half  angle)  should  constitute  a 
workable  design  for  atmospheric  aerosols. 
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